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Using an approach developed in the context of zero-temperature QCD to systematically 
sum higher order effects whose form is fixed by the renormalization group equation, we sum 
to all orders the leading log (LL) and next-to-leading log (NLL) contributions to the thermo- 
dynamic free energy in hot QCD. While the result varies considerably less with changes in 
the renormalization scale than does the purely perturbative result, a novel ambiguity arises 
which reflects the strong scheme dependence of thermal perturbation theory. 



I. INTRODUCTION AND SUMMARY 

The renormalization procedure in quantum field theory inevitably introduces a renormalization 
scale parameter into perturbative results; as /i^ is unphysical, its value is in principle arbitrary. 
Consequently, the calculation of any physical quantity would necessarily result in the dependence 
on /x^ disappearing. However, at any finite order of perturbation theory, residual dependence on 
jj? renders the result ambiguous, for changing the value of fj,'^ changes the predicted value of the 
physical quantity that has been computed [Q]. 

This problem has proved to be particularly acute in the calculation of the thermodynamic free 
energy in thermal field theory In quantum chromo dynamics (QCD) at temperatures much 
larger than the deconfinement temperature the free energy has been calculated ||3|, ^ ^ to order 

5 /2 

as , four terms in the perturbative series beyond the leading ideal-gas term. Whereas the first 

few approximations turn out to show little sign of convergence for any temperature of practical 

5/2 

interest, the result to order as happens to be centered about the results obtained in lattice gauge 
theory, but the dependence on the renormalization scale parameter fj? is so large that it has little 
predictive power. 

Recently, in Ref. Q in the context of standard model calculations, it was shown how the 
renormalization group (RG) equation can be used to sum in a systematic manner the leading log 
(LL), next-to-leading log (NLL), . . .effects given the perturbative results to one loop, two loop, 
. . . order. This so-called "renormalization group summation" (RGS) has been found to lead to 
a considerable reduction of the dependence on the parameter /x^ within a given renormalization 
scheme. 

In the case of thermal field theory this procedure requires generalization because the perturba- 
tive series also involves half- integer powers and logarithms of a^- This is carried out in Sec. and 
applied to the available three-loop result for the free energy in hot QCD in the modified minimal 
subtraction (MS) scheme. The RGE result which includes all LL and NLL effects turns out to 
be only weakly dependent on fj?. Unfortunately, as discussed in Sec. Ill, this does not increase 
substantially the predictive power of the three-loop calculation, because even within a fixed renor- 
malization scheme (here the MS scheme) there arises a new ambiguity in the form of the initial 
conditions for the differential equations, whose integration carries out the RG-summation. This 
ambiguity is somewhat larger than at zero temperature because the RG-summation of the thermal 



2 



perturbative series leads to two uncoupled sets of differential equations rather than one, reflecting 
the particular difficulties that thermal perturbation theory present. 

In the Appendix, we briefly discuss the difference between the strictly perturbative solution for 
the running coupling to two-loop order and the exact solution at two- loop order, which can be 
given in closed form in terms of Lambert's W function. 



II. RENORMALIZATION GROUP SUMMATION IN HOT QCD 



In QCD with nj flavors of quarks, the thermodynamic free energy at high temperature has been 
computed to be g, |, | 
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where p. is the renormalization scale parameter of the MS scheme and as{p) is the running coupling 
in this scheme whose form to three- loop order is given in though we shall restrict ourselves to 
its two-loop version in the following (see the Appendix for more discussion). 

Changing the renormalization scale parameter Jl in principle should not alter the value of F, 
as changes in the explicit fi are compensated for by changes of a<i(/i), and indeed, the result (|l|) is 

- 5/2 _ 

independent of /i to order as ■ However, numerically the dependence on ji is large, in fact larger 
than that of the result to order a], unless /i, which has to be of the same order as T to avoid large 
logarithms, is much larger than even the electroweak scale. 

Extrapolating from Eq. (|T|), we assume the complete all-order result of J- can be represented 
by a series of the form 



T/J^O = 1 + E [Rn{u)X- 
n=0 



where J^q is the ideal-gas value, x = ag (/x)/7r, u = xL,L = ln(/2V(2vrr)2) and 
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although only the coefficients with n < 1 are accessible by thermal perturbation theory. In fact, 
all of the perturbatively accessible cofficients with m = have been calculated already, with the 
exception of Cifi, which is forthcoming^. 

In Ref. it has been shown how to sum all RG-accessible logarithms when the lowest-order 
coefficients to the sums in Eq. (^, Anfi, Bnfl, and Cn,o, and the /3- function coefficients (see Eq. 
below) are known. In this paper we extend this to a perturbative series of the form (^). 



^ Y. Schroder, private communication 
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While the perturbative expression (|l|) is in powers of x^/^, successive RGE-perturbative expres- 
sions are given by 

^rL/^o = l + xRo{xL), (4) 
•^RGsZ-^o = l + xRo{xL) +x^/^So{xL), (5) 
•^RGsZ-^o = l + xRo{xL)+x^/'^So{xL) + x'^{Ri{xL) + To{xL)lnx), (6) 
•^RGsZ-^o = 1 + xRoixL) + x-^^^SoixL) + x'^{Ri{xL) + Tq{xL) Inx) + x^'^Si{xL) (7) 

which ah are perturbatively accessible in hot QCD, as we shall see. 

The explicit dependence of on and its implicit dependence through x{p?) are such that 



where 



P{x) = = {h2X^ + hx^ + b^x^ + ...)• 
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Substitution of (||) into (^) yields 
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Lowest order terms of the form x", x^^a and x" Inx give rise to differential equations for Ro{u), 
So{u) and Tq{u), respectively. The boundary conditions on these equations are the computed values 
of Aofi, Bo ft and Co,o respectively, which can be read off Eq. (|l|). This gives 



Ro{u) = Ao,ow-\ So{u) = BoflW-^/\ To{u) = Cq,^ 
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These functions incorporate the LL contributions to to all orders. 

To next order in the coupling x, we find differential equations for Si, and Ti (which rely on 
knowning the above solutions for Rq, Sq, and Tq). Solving these equations gives 
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which incorporate the NLL contributions to to all orders. 
Continuing in this way, we obtain 
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FIG. 1: Comparison of the renormalization-scale dependence of the RGS result for J- / J-q at T = 3Tc 
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for Uf = Q (sohd hne) and the perturbative result to order as (long-dashed line) when varying the 
renormalization scale /l around a central value of 2'kT by a factor . The two dots on the vertical axis give 
two recent lattice results from Refs. [|l2[ |l^. The short-dashed lines forming a big Z show the dependence 
of the RGS result on varying the arbitrary parameters k and by a factor of around 1. 
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these give the NNLL contributions to J^. 

The coefficients Aq^, Bofi, C'o.O) ^1,0, and Bi^ are determined in the MS scheme by Eq. (||) and 
this allows us to construct the RGS approximation J-j^Q^ given in Eq. (0). The next approximation 
would also involve Ti, R2, S2, and T2. For these, one would need results for Cifi and ^2,0) 
B2fl, C2,0) but only Ci^ is computable in perturbation theory because the order a^-contribution 
to J- is inherently nonperturbative ^. The remaining coefficients would have to be derived 
from a nonperturbative framework such as lattice gauge theory, or they might be estimated by 
nonperturbative resummation techniques such as Fade approximations ||lO| , O]. Given that, the 

(5) 

above formulae would allow the construction of the next RGS approximation ^^qj^- 



III. DISCUSSION 

In Fig. 1^ the renormalization scale dependence of the purely perturbative result ([l|) for T/J-q at 
T = 3Tc for 71/ = is displayed using a two-loop running coupling constant as(/i = e^2-KT) with 
A varying between —1 and 1 (long-dashed line). The two dots on the vertical axis give two recent 
lattice results from Refs. ||T^, and the perturbative result is seen to agree well with these for 
Ji = 2'kT (i.e. A = 0), but deviates strongly for different choices of the renormalization scale. 

The RGS result as obtained above is given by the solid line, and it shows a rather weak 
dependence on the renormalization scale. The systematic RG summation as described above is 
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thus able to absorb almost all of the scale dependence. 

However, this does not mean that the RGS approach predicts the perturbative result at 
p, = 2'kT. Rather, the latter has been used as the intial condition for the differential equations 
determining the functions i?o, 5*0, Tq, and Si. Organizing the RG summation in terms of a 
variable = / {k2t:T)'^] instead of the variable L introduced in Eq. (^) would have led to a 
different definition of the constants ^o,0; Bq^q, Co,o that provide the initial conditions. Moreover, 
the series involving half-integer powers of x in Eq. (|2|) leads to differential equations that decouple 
from those responsible for the series involving integer powers and logarithms of x. The former 
series could therefore have been derived by introducing a different variable Ly = / {v2'kT)'^] 
for Sn- 

In identifying the ambiguity in and Ljy (as parametrized by k and ly respectively) we are 
keeping the value of x unaltered — both Jl is held fixed and the form of x{jl) dictated by having 
chosen to work in the MS scheme is not changed. Rather, we note that the invariance of the 
perturbative result ([l|) under the change 

^1,0 + In -4^ ^ (^1,0 + ^1,1 In k) + In = a[^1 + Ai^iL^ (18) 

(with analogous equations for and c['^q) is lost when one deals with the RGS expressions 
Ri, Si and Ti. The boundary condition for the differential equations for Si and Ti, and 
the logarithm in the solution to these equations have a dependence on k (or v) that no longer 
automatically compensates as in the perturbative result. 

In Fig. H the short-dashed lines forming a big Z show the result of varying k and u around the 
value 1 by a factor with A between —1 and 1. The upper bar in the big Z is formed by = e^, 
K = e^, the diagonal hj u = k = e^, and the lower bar by = e~^, k = e^. Evidently, when k 
and ly are identified, the ambiguity is a little bit smaller than the one given by the renormalization 
scale dependence of the purely perturbative result, but varying independently of k leads to even 
larger variations, at their extremes even exceeding the ideal gas result for \ J^\. Unfortunately, when 
varying k and v, either together or independently, there is no saddle point for ^p^^ that would 
allow one to eliminate these ambiguities by a principle of minimal sensitivity [p. 

The renormalization group equation (^) states that the thermodynamic potential is independent 
of the renormalization scale p,. When one uses this equation to incorporate all the logarithmic 
contributions coming from higher order perturbation theory whose form is implied by this equation, 
one indeed finds the dependence of a perturbative approximation to T on jl being diminished, as 
expected. However, we have identified another source of ambiguity (characterized by k and i^) which 
leads to large variations of the perturbative RGS result for J^. This highlights the numerically 
strong scheme dependence of the perturbative result which still remains after having eliminated 
the strong dependence on the renormalization scale fl by the RGS method. (There are, of course, 
other ambiguities in the RGS result that could arise due to changing the renormalization scheme 
from MS; these we have not addressed.) 

As can be seen from Fig. ||, the extra dependence on the parameter u associated with the part 
of the perturbative series involving half-integer powers in ag is in fact the one which dominates 
the uncertainties of the perturbative result (the short-dashed lines are rather flat when only k is 
varied, but depend strongly on v). This part of the perturbative series is exclusively associated 
with "soft" collective phenomena such as screening and Landau damping and calls for a more 
complete treatment than conventional perturbation theory is able to achieve. Recent attempts 
in this direction that have been put forward include separate Fade approximations to soft and 
hard contributions WM, optimization of perturbation theory using the hard-thermal-loop effective 



action |15], and approximately self-consistent propagator resummation in the so-called <I>-derivable 
approach |16|. 



6 



Acknowledgments 

Financial support has been provided by NSERC. D. G. C. McK. would like to thank the Tech- 
nical University of Vienna where this work was performed. R. and D. MacKenzie provided helpful 
suggestions. Y. Schroder and G. Moore pointed out typographical errors in the first version of this 
paper. 

APPENDIX 



The beta function to two-loop order, 



2 dx dx 
d/u^ dL 



can be integrated in closed form 
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Hence, the solution to Eq. ( [A.l| ) can be written as a Lambert W function [^] with 
where the real branch W-i with W < —1 has to be taken. 
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FIG. 2: Comparison of the two-loop running coupling x — as/n resulting from the "exact" result ( |A.5| ) 
(solid line) with the strictly perturbative one from Eq. (A.6) (dashed line) for the case of QCD with n/ = 0. 
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On the other hand, the standard perturbative two-loop result is usually given as 

-1 / p lnL\ , 
x = -^ 1 + (A.6) 



which assumes small x and correspondingly large L. 

In Fig. |2| we compare the "exact" two- loop result ( |A.5D for x = as/vr (solid line) with the strictly 
perturbative one from Eq. ([A.6| ) (dashed line) for the case of QCD with Uf = (62 = —11/4, 
p = 51/22) as a function of L. The divergence at L = makes itself felt significantly earlier (as L 
approaches zero) in the standard perturbative result ( |A.6| ) than in the exact two-loop result ( A.5D . 



However, for L ^ 4, the difference between the two coupling is less than 1.5%. In the perturbative 
treatment of hot QCD, if one chooses a renormalization scale p = 2ttT this is indeed the case for 
all T > Tc ^ ^MS (™ pure glue QCD a typical value, which we have adopted in this paper, is 
Tc = 1.14Aj^); a noticeable difference thus arises only for smaller choices of p/T in combination 
with T sufficiently close to Tc. 
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